Reentrant Phase Behavior in Active Colloids with Attraction 
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Motivated by recent experiments, we study a system of self-propelled colloids that experience 
short-range attractive interactions and are confined to a surface. Using simulations we find that 
the phase behavior for such a system is reentrant as a function of activity: phase-separated states 
exist in both the low- and high-activity regimes, with a homogeneous active fluid in between. To 
understand the physical origins of reentrance, we develop a kinetic model for the system's steady- 
state dynamics whose solution captures the main features of the phase behavior. We also describe 
the varied kinetics of phase separation, which range from the familiar nucleation and growth of 
clusters to the complex coarsening of active particle gels. 



I. INTRODUCTION 

Active fluids composed of self-propelled units occur in 
nature on many scales, ranging from cytoskeletal fila- 
ments and bacterial suspensions to macroscopic entities 
such as insects, fish, and birds [1]. These systems exhibit 
strange and exciting phenomena such as dynamical self 
regulation [2], clustering [3-8], anomalous density fluc- 
tuations [9], and unusual rheological and phase behav- 
ior [10-13]. To elucidate the the microscopic origins of 
collective behavior, recent experiments have focused on 
developing simplified, controllable active fluids in nonliv- 
ing systems, using chemically propelled particles under- 
going self-diffusophoresis [14-17], Janus particles under- 
going thermophoresis [18, 19], and vibrated monolayers 
of granular particles [20-22]. 

Motivated by these experiments, we recently stud- 
ied [7] a minimal active fluid model: a system of self- 
propelled smooth spheres confined to two dimensions and 
interacting solely by excluded-volume repulsion. We dis- 
covered that activity in this system induces a continuous 
phase transition to a state in which a high density solid 
coexists with a low density fluid, complete with a bin- 
odal coexistence curve and critical point. The aggrega- 
tion of particles into dense clusters in the absence of any 
attractive interactions occurs as a result of self-trapping 
[4, 5, 7]. Separately, recent studies have examined sev- 
eral active colloidal systems which do possess attractive 
interparticle interactions [14, 23, 24]. The latter study 
demonstrates that activity suppresses phase separation, 
an effect which those authors postulate is generic. 

In this paper we investigate these opposing effects in a 
system of self-propelled particles endowed with both at- 
tractive interactions and nonequilibrium self-propulsion. 
We find that the phase diagram of such a system is reen- 
trant as a function of activity (Fig 1). This effect can 
be observed by traversing the phase diagram by increas- 
ing activity while holding the attraction strength fixed. 
At low activity the system is phase-separated due to at- 
tractive interparticle interactions, while increasing activ- 
ity to moderate levels suppresses phase separation and 
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FIG. 1. Top: Phase diagram illustrated by simulation snap- 
shots at time t — lOOOr with area fraction (f) — 0.4, as a 
function of interparticle attraction strength U and activity 
strength Pe. The colors are a guide to the eye distinguishing 
near-equilibrium gel states (green) from single-phase active 
fluids (red) and self-trapping-induced phase-separated states 
(blue). Bottom: Detail of reentrant phase behavior at U — 4 
as Pe is increased. At Pe = 4 (left), the system is nearly ther- 
mal and forms a kinetically arrested attractive gel. Increasing 
Pe to 20 (center) suppresses phase separation and produces a 
homogeneous fluid characterized by transient density fluctua- 
tions. Increasing Pe further to 100 (right) results in athermal 
phase separation induced by self-trapping. 



produces a homogeneous fluid. Increasing activity even 
further induces self-trapping which returns the system to 
a phase- separated state. We construct a simple kinetic 
model whose analytic solution captures the form of this 
unusual phase diagram and explains the mechanism by 
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which activity can both suppress and promote phase sep- 
aration in different regimes. We also describe the kinetics 
of phase separation, which differ significantly between 
the near-equilibrium and purely active phase-separated 
states. The behaviors we observe are robust to variations 
in parameter values, and thus could likely be observed in 
experimental systems of self-propelled, attractive colloids 
such as those studied in Refs. [14, 23, 24]. 



II. MODEL 



Our model is motivated by recently developed exper- 
imental systems of self-propelled colloids sedimented at 
an interface [14, 23], and consists of smooth spheres im- 
mersed in a solvent and confined to a two-dimensional 
plane [7]. Each particle is active, exerting a constant 
propulsion force on the surrounding fluid. Since the 
particles are smooth spheres and we neglect hydrody- 
namic coupling, they do not interchange angular momen- 
tum and thus there are no systematic torques. How- 
ever, the particles' self-propulsion directions undergo ro- 
tational diffusion; based on experimental observations 
[23], we confine the propulsion directions to be always 
parallel to the surface. For simplicity, interparticle in- 
teractions are modeled by the standard Lennard- Jones 

potential Flj =4e (^)^^ — (7)^ which provides hard- 
core repulsion as well as short-range attraction, with a 
the nominal particle diameter, and e the depth of the 
attractive well. 



The state of the system is represented by the posi- 
tions and self-propulsion directions {r^, of the par- 
ticles, and their evolution is governed by the coupled 
overdamped Langevin equations: 
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[i^Lj({rJ) + Fpi>,] 



(1) 
(2) 



Here Flj = — VVlj, Fp is the magnitude of the self- 
propulsion force, and i>i = (cos6>i, sin6>i). The Stokes 
drag coefficient 7 is related to the diffusion constant 
by the Einstein relation D = Dj^ is the ro- 

tational diffusion constant, which for a sphere in the 
low- Reynolds-number regime is I^r = The 77 are 

Gaussian white noise variables with {r]i{t)) = and 

{V^{tH{t'))=Si,5{t-t'). 

We non-dimensionalized the equations of motion using 
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a and /cbT as basic units of length and energy, and t = ^ 
as the unit of time. Our Brownian dynamics simulations 
employed the stochastic Runge-Kutta method [25] with 
an adaptive time step, with maximum value 2 x 10~^r. 
The potential Vi,j was cut off and shifted at r = 2.5cr. 



III. PHASE BEHAVIOR 



We parametrize the system by three dimensionless 
variables: the area fraction 0, the Peclet number Pe = 
'^pj, and the strength of attraction U = In or- 

der to limit our investigation to regions with nontrivial 
phase behavior, we fix the area fraction at = 0.4. At 
this density, a passive system (Pe = 0) is supercritical 
for U ^ 2.2, and phase- separated for stronger interac- 
tions [26]. For purely repulsive self-propelled particles, 
the system undergoes athermal phase separation as a re- 
sult of self-trapping for Pe ^ 85, and remains a homoge- 
neous fluid for smaller Pe [7]. 

To understand the phase behavior away from these lim- 
its, we performed simulations in a periodic box with side 
length L = 200 (with resulting particle count N = 20371) 
for a range of attraction strengths /7 G [1, 50] and propul- 
sion strengths Pe G [0,160]. Except where noted, each 
simulation was run until time lOOOr. Systems were ini- 
tialized with random particle positions and orientations 
except that (1) particles were not allowed to overlap and 
(2) each system initially contained a close-packed hexag- 
onal cluster comprised of 1000 particles to overcome any 
nucleation barriers. To quantify clustering, we consider 
two particles bonded if their centers are closer than a 
threshold, and identify clusters as bonded sets of more 
than 200 particles. The cluster fraction /c is then calcu- 
lated as the total number of particles in clusters divided 
by TV. 

The behavior of the system is illustrated in Fig. 1 by 
representative snapshots, and in Fig. 3 with a contour 
plot of /c. The most striking result is that the phase dia- 
gram is reentrant as a function of Pe. As shown in Fig. 1, 
low-Pe systems form kinetically arrested gels which grad- 
ually coarsen towards bulk phase separation. Increasing 
Pe to a moderate level destabilizes these aggregates and 
produces a homogeneous fluid, while increasing activ- 
ity beyond a second threshold accesses a high-Pe regime 
in which self-trapping [4, 5, 7] restores the system to a 
phase-separated state. 

As evident in Fig. 1, the widths of the low-Pe and high- 
Pe phase-separated regimes both increase with attraction 
strength [/, eliminating reentrance for U ^40. These 
trends can be schematically understood as follows. In the 
low-activity gel states, particles are reversibly bonded by 
energetic attraction. Particles thus arrested have random 
orientations, and so the mean effect of self-propulsion is 
to break bonds and pull aggregates apart. This opposes 
the influence of attraction, and so the width of the low-Pe 
gel region increases with U. By contrast, at high Pe we 
find that self-trapping is the primary driver of aggrega- 
tion. As shown in the next section, energetic attractions 
act cooperatively with self-trapping in this regime to en- 
able phase separation at lower Pe than would be possible 
with activity alone. 
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FIG. 2. Schematic representation of our kinetic model. The 
high- and low-density phases fill the left and right half-spaces. 
A particle on the cluster surface (red) escapes only if its di- 
rection points within the "escape cone" (green) defined by 
Fp(z> • n) > Fmax, while particles in the gas land on the sur- 
face at a rate proportional to the gas density and propulsion 
speed. 



IV. KINETIC MODEL 
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FIG. 3. Left: Fraction of particles in clusters fc measured 
from simulations as a function of Pe and U. Right: fc as 
predicted by our analytic theory (Eq. 3) which reproduces 
the major features of the phase diagram, including the gel 
region for Pe < [/, the self-trapping region for high Pe, and 
the low-/c fluid regime in between. The values of the fitting 
parameters are k, — 2 and /max — 1.7. 



To better understand the physical mechanisms under- 
lying the reentrant phase behavior, we develop a minimal 
kinetic model to describe the phase separated state. By 
analytically solving the model we obtain a form for fc 
which captures the major features of the phase behavior 
observed in our simulations. In the model we consider a 
single large close-packed cluster coexisting with a dilute 
gas which is assumed to be homogeneous and isotropic 
(Fig. 2). Particles in the cluster interior are assumed 
to be held stationary in cages formed by their neighbors, 
but their propulsion directions 6i continue to evolve dif- 
fusively. 

The rate per unit length at which gas-phase particles 
with a given orientation impinge on the flat surface of 



our model cluster is 



cos^, with the number 



density of the gas. Integrating over angles for particles 
traveling toward the cluster yields k[n = , the rate at 
which gas particles condense onto the cluster. 

Next we estimate the rate of evaporation. We note 
that a particle on the cluster surface remains bound so 
long as the component of its propulsion force along the 
outward normal, Fp{i> • n), is less than F^ax, the maxi- 
mum restoring force exerted by the binding potential. As 
shown in Fig. 2, this implies an "escape cone" in which 
the particle's director must point in order for it to escape. 

The critical angle is a = tt — cos~^ (^^^pi^)' ^^^^ /max 
the non-dimensionalized maximum restoring force scaled 
by the depth of the attractive well, which subsumes all 
relevant details of the binding force: /max = ^^W^k^' 

We now consider the angular probability distribu- 
tion of particles on the cluster surface P{6^t). If 
we take as our initial condition the steady-state dis- 
tribution P(^, 0) = P{0)^ then the rate per unit 
length at which particles leave the cluster is kout — 

The time evolution of P is 



A. 
at 



C^p{e,t)de 



t=0 



governed by diffusion, with absorbing boundaries at 
the edges of the cone: ^^^t^^ = oe^'*^ 



P{±a,t) 



0, 



with 



general solution P{0^t) = 
COS ^^^^ e~^'^ * , To simplify the analysis 

we note that higher-order terms decay rapidly in time, 
so the long-time behavior is dominated by the q = 1 
term. We therefore discard higher-order terms and ap- 
proximate the initial distribution as P{0) = ^ cos (from 
the distribution of incident particles) and solve to find 

^out — 4cra^ ' 

We further note from visual observation of our sim- 
ulations that the departure of a surface particle opens 
a hole through which subsurface particles may escape if 
their propulsion directions point outward. [7]. We des- 
ignate hz the total number of particles lost per surface 
particle that evaporates, and determine this value by fit- 
ting to our simulation data. The final expression for the 
evaporation rate is then /cout = 

Equating k^^ and /Cout yields a steady-state condition 
which can, by assuming a close-packed cluster, be rewrit- 
ten in terms of f^'. 



fc = 



(3) 



As shown in Fig. 3, this model reproduces the essen- 
tial features of our system, including active suppression 
of phase separation at low Pe, activity-induced phase sep- 
aration at high Pe, and a reentrant phase diagram. The 
model thus extends the analysis in Ref. [7] to describe 
the coupled effects of activity and energetic attraction. 
As noted in that reference, our model description of self- 
trapping can be considered a limiting case of the theory 
of Tailleur and Gates [4, 5] in which a self-propulsion ve- 
locity that decreases with density leads to an instability 
of the homogeneous initial state. 
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FIG. 4. Visual guide to phase- separation kinetics at fixed 
U = 30. To make mixing visible, particles are labeled in two 
colors based on their positions at t = 1. A passive system (top 
row) forms a space-spanning gel which gradually coarsens. 
The addition of activity (second row) greatly increases the 
rate at which the gel evolves. In both cases particles remain 
'local' and largely retain the same set of neighbors. When 
Pe exceeds U (third row), activity is strong enough to break 
the gel filaments and a fluid of large mobile clusters results. 
While the instantaneous configurations appear structurally 
similar to the gels above, these systems quickly become well- 
mixed due to splitting and merging of clusters (see also Fig. 
5). In the high-Pe limit (bottom row), self-trapping drives 
the emergence of a single well-mixed cluster surrounded by a 
dilute gas. 



V. PHASE SEPARATION KINETICS 

The kinetics of phase separation differ significantly be- 
tween the low-Pe gel and high-Pe self-trapping regions. 
In low-Pe systems, thermal influences dominate and the 
kinetics are those of a colloidal particle gel [27]. Since 
the area fraction in our simulations is high, the gels 
we observe appear nonfractal. Thermal agitation gradu- 
ally reorganizes the gel into increasingly dense structures 
[28, 29], leading toward a single compact cluster in the 
infinite-time limit. The presence of activity greatly ac- 
celerates the rate at which the gel evolves, as shown in 
Figs. 4 and 5. 

As Pe is increased beyond the threshold value Pe ^ [/, 
activity begins to overwhelm energetic attraction and the 
gel is ripped apart. This arrests the compaction, result- 
ing in a plateau in the system's total potential energy 
(Fig. 5). Just above this transition, the system resem- 
bles a fluid of large mobile clusters which rapidly split, 
translate, and merge (Fig. 4). This phase is clearly iden- 
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FIG. 5. Quantitative measurements of states with U — 30. 
Left: Mean interaction energy per particle {E represents the 
total system potential energy). Low-Pe gels are strongly ar- 
rested and do not approach bulk phase-separation on our sim- 
ulation timescales; however, higher-Pe systems evolve much 
faster and nearly reach the bulk limit. When Pe > U (dashed 
lines), coarsening is arrested as the gel breaks into a fluid of 
mobile clusters with a Pe-dependent characteristic size, lead- 
ing to a plateau in the system energy. Right: Discrimination 
of gel from fluid states by mean-square displacement. Dotted 
gray lines have slope 1 and highlight the distinction between 
subdiffusive gels (Pe < U), and superdiffusive active fluids. 



tified by mean-square displacement measurements (Fig. 
5) which are superdiffusive for this active fluid, unlike 
the subdiffusive behavior found in gels. As Pe is further 
increased, the characteristic mobile cluster size decreases 
until the appearance of an ordinary single-phase active 
fluid is recovered. 

Additional increase of Pe will eventually cross a sec- 
ond threshold into a phase-separated regime whose be- 
havior is dominated by self-trapping. As reported previ- 
ously [7], these systems undergo nucleation, growth, and 
coarsening stages in a manner familiar from the kinet- 
ics of quenched fluid systems, albeit with the unfamiliar 
control parameter Pe instead of temperature. 

These three regimes are characterized by dramatically 
different particle dynamics. To illustrate these behaviors 
and their effect on particle reorganization timescales, in 
Fig. 4 we present snapshots from simulations in which ini- 
tial particle positions are identified by color. We see that 
when attraction is dominant {U > Pe), each particle's 
set of neighbors remains nearly static over the timescales 
simulated, indicative of long relaxation timescales due 
to kinetic arrest. In contrast, when activity dominates 
(Pe > U) particles rapidly exchange neighbors and the 
system becomes well-mixed. Importantly, note that par- 
ticle dynamics cannot be directly inferred from the in- 
stantaneous spatial structures in the system. For exam- 
ple, while systems with low activity (Pe < [/, Fig. 4 
second row) and moderate activity (Pe ^ [/, Fig. 4 third 
row) appear structurally similar, the rate of particle mix- 
ing differs by orders of magnitude. 
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VI. CONCLUSIONS 



Activity can both suppress and induce phase separa- 
tion, and we have shown that these opposing effects can 
coexist in the same simple system. The resulting coun- 
terpoint produces a reentrant phase diagram in which 
two distinct types of phase separation exist, separated 
by a homogeneous fluid regime. This surprising result 
makes it possible to use two experimentally accessible 
control parameters (Pe and U) in concert to tune the 
phase behavior of active suspensions. This control is es- 
pecially valuable because attractive interparticle interac- 
tions are common in experimental active systems, being 



either intrinsic [14, 23] or imposed, such as by the addi- 
tion of depletants [24] . An understanding of the complex 
phase behavior accessible to these systems is a critical 
stepping stone toward designing smart active materials 
whose phases and structural properties can dynamically 
respond to conditions around them. 
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